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Abstract. Turbulent coagulation in protoplanetary disks is known to operate on timescale far shorter than the 
lifetime of the disk. In the absence of mechanisms that replenish the small dust grain population, protoplanetary 
disks would rapidly lose their continuum opacity-bearing dust. This is inconsistent with infrared observations 
of disks around T Tauri stars and Herbig Ae/Be stars, which are usually optically thick at visual wavelengths 
and show signatures of small (a < 3/im) grains. A plausible replenishing mechanism of small grains is collisional 
fragmentation or erosion of large dust aggregates, which model calculations predict to play an important role in 
protoplanetary disks. If optically thick disks are to be seen as proof for ongoing fragmentation or erosion, then 
alternative explanations for the existence of optically thick disks must be studied carefully. In tis study we explore 
two scenarios. First, we study the effect of residual, low-level infall of matter onto the disk surface. We find that 
infall rates as low as 10 _11 MQ/yr can, in principle, replenish the small grain population to a level that keeps the 
disk marginally optically thick. However, it remains to be seen if the assumption of such inflow is realistic for 
star+disk systems at the age of several Myrs, at which winds and jets are expected to have removed any residual 
envelope. The effectiveness of even a low level infall can be understood by the strongly non-linear behavior of the 
coagulation equation: a high, fine-grain, dust density at any given time leads to very, effective removal of these 
small grains, while a low fine-grain density lasts for a far longer time. We then consider a second scenario in which, 
during the buildup phase of the disk, an intermediate fine-grain dust abundance is generated that is sufficiently 
low to ensure longevity of the state yet sufficiently high for the disk to remain optically thick. While our models 
confirm that such an "initial condition" can be constructed, we argue that these special initial conditions cannot 
be achieved during the disk build-up phase. In summary, fragmentation or erosion still appear to be the most 
promising processes to explain the abundant presence of small grains in old disks. 
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1. Introduction 

The first stages of planet formation are now directly 
observable with infrared observations of protoplanetary 
disks. Most studies focus on the presence of slightly 
grown dust gr ains as evidence o f the o n set of this growth 
process Ce.g. Ivan Boekel et all 2003; P rzygodda et al 



20031: Ivan Boekel et alj . 12004 iKessler-Silacci etal.l . l2006t ). 
However, a much undervalued fact is that even old 
(~3 Myr or even older) disks still contain significant 
populations of small (< 3/im) grains as indicated by 
the pres ence of stron g infr ared solid state emission 
features dMalfait et al. , 1998) and large optical depth 



band, which causes the disk to absorb or scatter all 
incident stellar flux Models of dust coagulation with- 
out fragmentation predict consistently severe depletion 
in these small grains on timescale far s maller than the 
disk l ife time, as seen in calcul a tions by IWeidenschillinel 
(Il984l). iDullemond fc Dominikl (j2005l hereafter DD05), 
and lTanaka et all (|2005l hereafter TIH05). While the as- 
sumption of a perfectly laminar disk may allow a small 
population of micron-size grains to survive (TIH05) , even 
very weak turbulence (at a level of a = 10~ 6 ) would en- 
hance the coagulation process sufficiently for the disk to 
become optically thin. In the above models, differential ra- 
dial drift was not included: at low levels of turbulence, it is 
the dominant mechanism controlling the relative velocities 
of grains of various sizes ([Whipple! Il972t IWeidenschillinel . 



19771 ll997tlBrauer et all 120081 ), which implies that, with 



( Bouwman et aLfliool) . By optically thick, we refer out fragmentation, growth is enhanced by this mechanism. 



to a significant optical depth as measured in the V 



2 



Dominik, Dullemond: Coagulation of small grains in disks 



Models of dust coagulation without a treatment 
of collisional f r agmentatio n are of course unrealistic. 
Weidenschillina (|l984ll997l ) included fragmentation in his 
model. It turns out that fragmentation is a universal out- 
come. Models of particle motion in laminar and turbulent 
accretion disks predict collision velocities in excess of the 
critical value for break-up (Weidenschilling 1977; Volk et 
al. 1984; Blum & Wurm 2000; for a review, see Dominik 



et al. (2007)1. iRrauer et all (|2008l ) demonstrated that col- 
lisional fragmentation is a major barrier to the growth 
of dust aggregates beyond centimeter or decimeter sizes. 
However, it is still unknown if fragmentation can produce 
a sufficient amount of small grain dust to produce the 
disk opacity inferred from observations. This question is 
difficult to answer, because the amount of small grains 
produced in collisions is strongly dependent on collision 
parameters (Paszun and Dominik, submitted to A&A). It 
is therefore important to know if the opaqueness of proto- 
planetary disks is an indication of small grain production 
by fragmentation and erosion (as we argued in DD05), or if 
other mechanisms could retain the small grain population 
in disks, such as the thermal breakup of icy aggregates 
(Tanaka et al. 2005, poster at PPV), which can operate in 
a limited radial range. 

In the present paper, we investigate if a continuous 
replenishment of small grains by residual infall of mat- 
ter from the circumstellar envelope onto the disk's surface 
could sustain the small dust grain population. This in- 
fall might be residual infall from the paren t al clo ud, or 
Bondi-Hoyle like accretion (jBondi fc Hovlel . Il944h by a 
protoplanetary disk moving through a molecular cloud. 
We assume that the infall originates in the parental cloud, 
with zero velocity at infinity. As long as the motion of the 
star relative to the cloud is small compared with Keplerian 
velocities in the inner regions of the disk, these cases are 
simila r. The idea of e xtern al replenishment was first stud- 
ied bv lMizuno et al. ( 1988), who were interested primarily 
in convection-driven turbulence (which requires optically- 
thick disks) during the early stages of the disk's lifetime. 
Using a simple one-zone model of a disk, they indeed 
showed that a mass infall rate of 10~ 6 and 10- 7 M Q /yr 
could help to keep the optical depth of the disk high. 

One can easily show that an envelope with an infall 
rate of 10~ 6 M Q /yr would be optically thick and com- 
pletely obscure the star+disk system, while 10 _7 M Q /yr 
would still be marginally optically thick. The typical 
star+disk systems studied in the context of grain growth 
and planet formation are 'revealed sources' (class II in 
the Lada classification), meaning that they have envelopes 
that are optically thin. If this process is to work at all in 
these systems, it must work at relatively low infall rates, 
which we consider here. 

In Sect. [3] we first study if the replenishment of the 
disk with small dust grains by an infalling tenuous en- 
velope with low infall rate would be sufficient to ensure 
that the disk remained sufficiently optically thick to be 
consistent with observations. We employ the 1-D vertical 
slice models for coagulation described by DD05, but with 



additional continuous infall of small grained dust at the 
top. We study the total vertical optical depth, the optical 
depth in grains smaller than 3 /im (which are the silicate- 
feature-producing grains) , and the optical depth in grains 
larger than 3 /jm. 

In Sect. [J] we consider the effects of a low initial abun- 
dance of small grains on the removal timescale of small 
grains. This involves models without external inflow, but 
with different dust-to-gas ratios in the initial setup. 

We discuss our results and present our conclusions in 
Sect. [S] 

2. The model 

2.1. Disk and coagulation model 

We model the dust coagulation in a disk around a pre- 
main-sequence star. We begin the simulation at t = 0, 
which represents the time immediately after the main 
build-up phase of the star+disk system. This is the mo- 
ment when the main infall onto the system, which is of 
the order of 10~ 5 Mo/yr during the build-up phase, either 
ceases entirely or reduces to far lower values. It is this 
much lower 'residual infall' that we focus on in Sect. [3l 

As our base star+disk model, we consider a 0.05M Q 
disk around a 0.5M Q star. The star is assumed to have an 
effective temperature of 4000 K and a radius of 2.5R Q . For 
the surface density of the disk, we assume a power law S = 
So( r / AU) -1 (see e.g. Dullemond et al 2007 for a review) 
where r is the distance from the star. The disk radius is as- 
sumed to be 200 AU, and M disk = 0.1M* = 0.05M Q , such 
that the surface density at 1 AU becomes So = 354g/cm 2 
(in gas). The vertical gas distribution is assumed to be in 
hydrostatic equilibrium, with a gas temperature given by a 
simple irradiation recipe T gas (r) = (l/V2)4> 1 / 4 T*^R±/r 1 
where <fi is the irradiation angle, which we assume to be 
0.05. The dust and gas are initially well-mixed, with a 
dust-to-gas ratio by mass of 0.01. The dust grain size dis- 
tribution in the initial disk model and in the infalling 
flow of matter is assumed to be an MRN distribution 



(Mat his et all 119771 ) with particles distributed between 



0.1 and 0.5^m, following a power-law size distribution 



/(a) 



tx a 



-7/2 



The dust grains are assumed to be as- 



tronomical silicate. The opacities were c omputed using 
Mie t heory, with optical properties from lLaor &: Drainel 
(Il993l) . 

The coagulation model is that of DD05, which we used 
unmodified for these simulations. The disk's gas density 
and temperature were assumed to remain constant in spite 
of the vastly changing dust properties with time. We in- 
cluded coagulation by Brownian motion, vertical settling, 
and turbulent motions. The calculations do not include 
coagulation driven by radial drift in the disk. Radial drift 
is known to be im porta nt for the growth of large particles 
(IWeidenschillind . Il980l Il984l ). *n the current paper, we 



are interested in the small grain component, for which the 
differences are expected to be minor. Furthermore, it al- 
lows us to compare directly with the results by DD05, who 
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Table 1. Values for various quantities at different distances r from the star, for the model with trickling rate 
10 -9 MQ/yr. The quantities are: surface density £, trickling flux /trickle, impact velocity of the infall vg, midplane tem- 
perature T m id, pressure scale height H p , accretion shock height H B h, stopping time for 0.1/xm grains i s top, deceleration 
length for 0.1/zm grains in units of the shock height. 
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emphasized a need for some replenishment mechanism of 
small grains. In a model with strong turbulent mixing, the 
coagulation occurs particularly rapidly as material from 
the disk upper layers is mixed more deeply inside the disk, 
where the densities are higher and coagulation timescales 
are far shorter. 

2.2. Infalling envelope 

For the computations in this paper, we assume that the 
infalling matter is the residual influx from the collapse of 
a rotating cloud. The density and velocity distribut ion for 
this ca se were described with analytical models by lUlrich 
( 19761 ). represent ing a simplified ana lysis of a the more 
general study by iTerebev et al.1 ( 19841 ). The main param- 
eters of this model are the infall rate M and the centrifugal 
radius r contr , which we fixed to be 200 AU. The gas density 
at a distance r and polar angle fi = cos 9 is given by: 



p(r,6) = 



M 



i + ii 



4WGMr 3 \ Mo 
where /zq is a solution of the equation: 



Mo r 



— i-f)=l-M3< • 

^centr \ MO , 



The r- and 0-velocity components are: 

1/2 / x 1/2 

M 

Mo, 




(1) 
(2) 

(3) 
(4) 



Multiplying v s (9 = tt/2) by the density p(9 = 0/2) pro- 
vides the mass flux onto the disk surface: 



M 



trickle 



87rr/i r ce ntr 



(5) 



where = \/l - r/r ccntr . 

Some example values for the trickling inflow onto the 
disk surface are given in Table [1] 



2.3. Location of the accretion shock and dust injection 

As matter falls onto the disk, an accretion shock must form 
to decelerate t he gas. The physics of th e accretion shock is 
complex (e.g. iRuzmaikina fc id . Il994l ). because the tem- 
perature behind the accretion shock has to be computed. 
In the current study we assume that the cooling behind 
the shock is fast, i.e. that the post-shock gas quickly cools 
to the local temperature of the disk. We may then compute 
the location of the shock by equating the ram pressure in 
the infalling material to the local hydrostatic pressure in 
the disk and solve for the height above the disk. The ram 
pressure from the infalling material is given by 



Mr 



p(r,0H 2 (r) 



(6) 



The hydrostatic pressure in the disk pdisk is given as a 
function of height z in the disk by 



Mdisk = 7PdiskC s 



H„V2tt 



exp 



(7) 



where H p = c s /f2fc is the pressure scale height in the disk, 
c s = y/kT//j,m p is the sound velocity, jj, is the mean molec- 
ular weight, m p the proton mass, f2& = \/ GM*/r 3 is the 
local Kepler frequency, and 7 is the adiabatic index. 



Zsh 



\ 



-2H p \n 



(8) 



We find that this is typically 5 to 6 pressure scales above 
the disk (for detailed values, see Table [T]). 

When the infalling gas motion has stopped, the dust 
particles are injected into the disk with the free-fall speed 
at the considered location. Dust grains are decelerated 
over a certain distance, which we compute numerically. 
The results are reported in Table [T] and indicate that 
the stopping time is approximately a few hours, and the 
stopping length is far smaller than the disk scale height. 
Another issue is the temperature of the grains during de- 
celeration. The grains heat up due to frictional coupling 
with the gas. We can estimate the heating with simple ar- 
gument. The kinetic energy of the incoming dust grain is 
given by 0.5TO gr t>g. Following the deceleration curve of the 
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grain, we make the extreme assumption that the entire ki- 
netic energy is transformed into heating grains (in reality, 
a large fraction will heat the surrounding gas). Equating 
this heating rate to a radiative cooling rate, we found that 
the grains were heated up by less than a few hundred de- 
grees at 1AU, and that the heating was negligible at 10 
and 100 AU. The grain should therefore be able to reac h 
the disk largely unchanged fe.g. iRuzmaikina fc id , 119941 ). 

We can therefore justify our simple treatment of dust 
injection. Dust particles are added to our uppermost grid 
cell in the local disk slice with an initial velocity of 0. From 
that location, particles will settle down and become part 
of the coagulation process. 

3. Results I: Effect of residual infall 

We completed the coagulation simulations at three differ- 
ent radii in the disk: 1 AU, 10 AU and 100 AU. We as- 
sumed three different infall rates (0, 10 -11 , 10 _7 )M Q /yr, 
which translate into local infall rates in g cm~ 2 s _1 using 
Eq. ((5]). We also assumed two different strengths of turbu- 
lent mixing: a weakly turbulent disk (a = 10~ 6 , as a nearly 
laminar case) and a strongly turbulent disk (a = 10~ 2 , as 
typically found by MM instability calculations, and also 
needed to explain the high accretion rates during the early 
stages of disk evolution). Here, a is the Shakura-Sunyaev 
parameter for viscosity, which can be used to parametrize 
the strength of the turbulence, assuming a Stokes number 
of unity. We do not include fully laminar conditions, since 
we believe that these are unrealistic. We allow each model 
to run for 1 Myr, in a few time-consuming exceptions we 
stopped earlier if a steady state had already been reached. 
This provides a 3x4x3 matrix of models. 

In Fig. [2l the resulting time-dependent optical ver- 
tical depths at A = 0.55 /im are shown for the case of 
a = 10~ 6 . The top row shows the results without resid- 
ual infall, which can be most closely compared with the 
laminar models of DD05 and TIH05. It should be noted, 
however, that the small, but non-zero, turbulence in the 
present models has non-negligible effect on the results, and 
the comparison is therefore only partly valid. In these top 
three models, the coagulation is first driven by Brownian 
motion, which is responsible for the weak reduction in op- 
tical depth between t —0 and t —300 yrs for the R=1AU 
case. The optical depth then suddenly begins to drop, as 
sedimentation-driven coagulation reaches its peak inten- 
sity and efficiently removes small grains from the disk. 
The slow further decline in optical depth is then caused by 
turbulence-driven coagulation, which slowly removes ap- 
proximately 0.1 percent fine dust that escaped the sweep- 
up by settling grains. The final optical depth at 1 Myr 
turns is 0.2 at 1 AU, 2 at 10 AU, and 12 at 100 AU. While 
the outer regions of the disk are still optically thick, the 
inner regions have become optically thin. This was found 
by TIH05 to occur in their laminar models and they pro- 
posed that this might be a possible me chanism for explain- 
ing 'inner holes' i n certain disks (e.g . iForrest et all 120041 
Bouwman etaD . l2003t ICalvet et all |2002j). The second 



row of Fig. Q] shows the same models, but with a continu- 
ous trickling infall at a (global) rate of M — lO _11 M /yr. 
We see that at 1 AU this causes the optical depth to be- 
come constant at a value of 2. The outer disk is hardly 
affected. At an infall rate of M = 10 -7 M Q /yr, lower row, 
the optical depth at 1 AU becomes constant at 80 and 
at 10 AU and 100 AU at 70. These values are sufficiently 
high for the disk to be considered optically thick disk. 

Figure[5]shows the same models as Fig.[Tjbut now with 
a = 10~ 2 . Again, the top row shows the models without 
infall. As for the low-a case, Brownian motion produces an 
initially slow decline in the optical depth up to the onset 
of sedimentation-driven coagulation. However, in contrast 
to the low-a case, the subsequent decline in optical depth 
is extremely fast, which differs from our results in Fig. [TJ 
Turbulence-driven coagulation is far stronger than in the 
a = 10~ 6 case, and becomes highly effective in removing 
small grains. In these simulations we emphasize that the 
optical depth at late times is due primarily to particles 
larger than 3 /im, an effect that might be artificial and 
in part caused by the neglect of grain drift, unlike the 
low-a models. In the second row, we see that at 1 AU 
and 10 AU the trickling infall of fresh matter leads to a 
slow recovery in the optical depth, after its strong previous 
decline. The optical depth is now primarily controlled by 
small particles. At 1 AU, the optical depth levels off at 
around 1, and at 10 AU becomes 0.1, which is still optically 
thin, and lower than at 1 AU. At this radius, the optical 
depth is still increasing after 1 Myr. At 100 AU the infall 
has not produced any noticeable effect. In the last row, 
we can see that following the recuperation of the optical 
depth after its reduction it appears to saturate at a value 
of 150 at 1 AU as well as at 10 AU. The system reaches an 
equilibrium state, in which the rate of fine dust removal 
appears to equal the rate of replenishment by fresh fine 
dust. The equilibrium state is still imperfect, since the 
large grain population continues to oscillate, as can be 
seen by the dashed lines. However, since these large grains 
(> 3/im) represent only a small contribution to the optical 
depth, this is of no relevance to the optical depth issue. 
Even at 100 AU, we now see that the replenishment by 
new grains is effective and increases the optical depth, 
but an equilibrium state is not yet reached. 

Judging from these results, we see that, in principle, 
even a small residual infall rate onto the disk can explain 
the observed optical thickness of T Tauri star disks, even if 
fragmentation is not active and therefore cannot replenish 
the population of small grains. Whether infall scenarios 
with these rates in the late evolutionary stages are realistic 
is, however, still a major issue of concern, which we discuss 
in Sect. [5] 

4. Results II: Effect of reduced initial fine-dust 
abundance 

In the models of coagulation described here, it is seen that 
the dust in the disk coagulates initially extremely rapidly, 
reducing the optical depth to almost zero apart from that 
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r = 1 AU r = 10 AU r = 1 00 AU 




10° 10' 10 2 10 3 10 4 10 5 10 e 10 1 10 2 10 3 10 4 10 5 10 6 10 1 10 2 10 3 10 4 10 5 10 6 
Time [yrs] Time [yrs] Time [yrs] 

Fig. 1. The vertical V-band optical depth of the disk with turbulent parameter a = 10~ 2 at 1 AU, 10 AU and 100 
AU (from left to right) and infall rates 0, 10~ Mo/yr, lO -7 M0/yr (from top to bottom). The dotted line is the 
contribution to the optical depth by grains with radius < 3/im, the dashed line is the contribution by grains > 3/im. 



due to a trickling accretion of new material. It is striking, 
however, that even a very low residual infall rate of, say, 
10 _11 MQ/yr can have this effect. At first sight this ap- 
pears difficult to understand, because at this trickling rate 
it takes 5 Gyr to replenish fully a disk of 0.05M Q , whereas 
the removal of small grains occurs within 10 4 years. One 
would expect the replenishment to be unable to compete 
with the coagulation. The speed of coagulation is however 
inversely proportional to the density of dust. Therefore, 
the time to grow to sizes that no longer contribute to the 
optical depth will be inversely proportional to the density 
of fine-grain material. If we begin with more fine-grained 
dust, we therefore end with less. 

This suggests another way to retain the value the op- 
tical depth in a disk for a long time. Our simulations have 
consistently started from an already constructed disk, 
i.e. after the early build-up phase of the disk. Without 
justification, we have assumed that at this point in time 
the dust population is interstellar in both size distribu- 
tion and abundance. However, coag ulation is also active 
during the build-up pha se of the disk (jSuttner et all fl999: 
Suttner fc Yorkd . 120011 ). and it is therefore to be expected 



that by the time the main infall phase is over (and our 
simulation begins) the abundance of fine-grain dust has 
already decreased. 

To investigate the significance of this reduction, we 
could apply the models of Sect. [3] to the build-up phase, 
and observe the equilibrium state that is reached. We sim- 
ply need to adopt a far higher infall rate, i.e. a rate at 
which a star can be built in about 0.1 Myr. In our case 
this corresponds to an infall rate of about lO -5 M0/yr. In 
Fig. [31 we see that, at these high infall rates, the equilib- 
rium is reached at a very early stage, justifying the use of 
this state as an initial condition. In this particular case the 
figure shows that the optical depth is reduced by about 2 
orders of magnitude compared to the "interstellar" dust 
abundance case. As an initial condition of our Myr-long 
simulation, it would therefore be justifiable to begin with 
a fine-dust abundance that is a factor of 100 below the 
interstellar medium value. However, we have not modeled 
this buildup phase consistently, and more sophisticated 
physical ingredients, such as pre-collapse coagulation, disk 
accretion, and radial drift of larger grains, would need to 
be included, which is well beyond the scope of this paper. 
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10° 10' 10 2 10 3 10 4 10 5 10 e 10 1 10 2 10 3 10 4 10 5 10 6 10 1 10 2 10 3 10 4 
Time [yrs] Time [yrs] Time [yrs] 

Fig. 2. Like fig. [2] but with with turbulent parameter a = 1CP 6 . 
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We therefore study the effect of a reduced initial fine-dust 
abundance on the Myr-long simulations, by using the ini- 
tial abundance as a free parameter. This is implemented 
by completing models without any residual infall, but with 
decreased initial dust-to-gas ratio. We ignore any popula- 
tion of very large grains that may have been produced 
during the build-up phase. 

The results displayed in Fig. |4] indicate that the initial 
optical depth is proportional to the initial dust-to-gas ra- 
tio (which is obvious) and that this optical depth can be 
sustained for a time that is inversely proportional to the 
initial dust-to-gas ratio. For example, with an interstellar 
dust-to-gas ratio (i.e. 10~ 2 ) the optical depth is almost 
10 5 at first, but a steep decrease occurs after 100 years. If 
we begin the simulation, however, with a dust-to-gas ratio 
of 10~ 7 (i.e. a factor of 10 5 lower than interstellar), then 
the initial optical depth is close to 1 and can be sustained 
for about 1 Myr. It therefore appears that the strong re- 
duction of optical depth may occur as late as 1 Myr, but 
only if start with an optical depth that is already close to 
1. 

This implies that there exists a small region of parame- 
ter space in which the disk can remain marginally optically 
thick (ry ~ 1) for 1 Myr, without fragmentation or any 



other means of replenishment of small grains. However, 
many disks live even longer, and there is a large vol- 
ume of observational evidence suggesting that these disks 
have higher optical depths (e.g. Forrest et al. (2004) , and 
Dullemond and Dominik (2004) ; for a review of fitting 
disk models to observations, see Watson et al, (2007)). 

It will also be difficult to produce this specific initial 
state, for the following reason: low optical depth equilib- 
rium states require a long time to be established. The 
equilibrium state reached during the disk build-up phase 
is defined to be the state in which the timescale of re- 
plenishment of fine dust equals the timescale of removal 
of fine dust by coagulation. Therefore, the build-up time 
of the state with replenishment equals the lifetime of the 
state without replenishment. For example, during the disk 
build-up phase (clearly at a low infall rate) we would re- 
quire 1 Myr to achieve an equilibrium state that could, 
after the end of infall, be sustained for another 1 Myr 
as an optically thick disk. It is believed, however, that 
the build-up phase of a star+disk system (i.e. the class 
O and class I phases in the Lada/ Andre scheme) has a 
duration of les s than 1 Myr ( Andre fc Montmerlel Il994t 
Barsonv , 1994}) . implying that this phase is significantly 
shorter than the typical disk lifetime. The conclusion is 
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Fig. 3. The time-dependent optical depth for the case 
when we apply our model to the build-up phase of the 
disk, with an infall rate of 10 _5 MQ/yr. Since this sim- 
ulation is merely meant as an 'initial condition' for the 
Myr-long coagulation process studied in this paper, we 
run it only until an equilibrium state has set in. As can 
be seen in this figure, this happens after about 10 4 years. 



5. Discussion 

5.1. Consequences for dust fragmentation 

In Sect. 21 we have shown that a low value of fine-grained 
dust density in the disk can in principle persist for a long 
time, but that there is only a narrow range in parame- 
ter space for which the disk optical depth remains above 
unity for a few million years. It also takes a long time 
to reach this long-lasting state. This result is fully inde- 
pendent of the mechanism that produces the small grains 
in the disk. We can also apply it to small grains pro- 
duced by aggregate fragmentation and erosion (instead of 
trickling infall). If, after studying possible alternatives in 
this paper, we return to fragmentation and erosion as the 
most likely sources of small grain replenishment in disks, 
then we must conclude that these must be ongoing pro- 
cesses lasting throughout the few-million-year lifetime of 
the disk. If this replenishment process were to switch off 
for some reason, the disk would quickly become optically 
thin. 

5.2. Limits on infall onto T Tauri star disks 




Time [year] 

Fig. 4. The time-dependent optical depth for the for disks 
with varying initial dust-to-gas ratio. The curves are la- 
beled with log/dg. This figure clearly demonstrates the 
"The higher you start, the lower you end" effect caused 
by the non-linear behavior of the coagulation equation. 



How realistic is the assumption of low-level infall onto 
the disk? There are several ways to place constraints us- 
ing both observational and theoretical arguments. In the 
following, we show that acceptable infall rates are con- 
strained to be in the narrow range between 10 -7 and 
l(r 8 M Q /yr. 

5.2.1. Optical depth of the envelope 

We emphasize hat we consider sources with only lim- 
ited extinction toward the star. Many of the sources 
with optically thick disks have circumstellar reddening 
Ay lower than (sometime s far lower than) two (e.g. 



van den Ancker et al.l . 119971 ) . The largest contribution to 
the optical depth of the infalling envelope originates in the 
region around r ~ r ccntI . Using the density distribution 
from the Ulrich infall model, we arrive at an optical depth 
in the visual of 1.88 for an infall rate of 4 x lO _7 M0/yr. 
Therefore, the maximum infall rate that is consistent with 
the class II status of these sources is about this value. The 
optical depth at different infall rates obviously scales lin- 
early with the infall rate. 



5.2.2. Stellar winds 



that the dust cannot reach a steady state that lasts for 
a few Myrs. Any steady state reached during this time is 
certain to decay shortly after the main infall phase. It is 
also possible that no steady state of the dust population 
is reached at all during this disk build-up phase. Both sit- 
uations are not promising for keeping the optical depth 
high over a few Myrs time. 



Young stars are also known to possess stellar winds, which 
can remove residu al circumstellar enve lopes. Using in- 
frared and UV lines, Dupree et alJ ( 2005h argued that TW 
Hya accelerates a hot wind with a velocity of 400km/s and 
a mass loss rate of about lO _11 M /yr, which should be 
able to quench any infall with rates below 10 -9 M Q /yr. 
The evidence for this wind w as however questioned by 
Johns-Krull fc Herczed ( 2007 ). If these wind properties 
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were typical for T Tauri stars, 10 9 M Q /yr would be a 
lower limit for residual infall to reach the disk. 



5.2.3. Radiation pressure and dust/gas separation 

Another obstacle to any residual infall is radiation pres- 
sure. At low infall rates, the envelope is optically thin and 
the radiation reaches all the dust grains in the infalling 
gas. The problem of radiative blowout of grains from 
main-sequence stars was studied extensively in the con- 
text of the Veg a phenomenon and for meteorites in the 
Solar System dBurns et all Il979[ bruen fc Zookl . Il980t 
Artvmowicz , 19881 : Krivova et all 200C ) . For the luminos- 
ity of the Sun, a small range of dust particle sizes around 
0.1/Ltm have the correct prop erties for rad i ation pressure to 
exceed gravitational forces ( Burns et al. . 19791 ). For more 



luminous sources such as A stars or even their still more lu- 
minous Herbig Ae precursors, the blowout limit is pushed 
to 10/im and beyond, such that grains in a size range 
from 0.01 to lOjum should be blown away from the star 
( Artymowiczl . Il988l ). The one difference in the scenario 
considered here is of course that the dust grains are em- 
bedded in an infalling gas envelope, which is not affected 
by any significant radiation pressure and could possibly 
drag the dust particles along. 

The drift velocity caused by the radiation pressure 
on dust grains must then exceed the infall velocity of 
the material. We can place a limit on the mass accre- 
tion rates necessary for dragging the grains along in the 
following way. In a spherical free-fall accretion scenario 
with a constant infall rate M, the local velocity is given 
by v(r) = —\/2GM±/r and the gas density is given by 
p = M / {Airr 2 \v\). We assume that a dust grain is embed- 
ded in that infall. The stellar radiation in incident on the 
grain with a flux L+/ (Airr 2 ). The radiation force on the 
grain is then 

T-i 2 



47Tr 2 ( 



(9) 



where a is the grain radius and c is the speed of light. 
The efficiency for pressure transmission to the grain is 
Qpr = Qabs + Qscat(cosa) is , with Q abs being the ab- 
sorption efficiency, Q S cat the scattering efficiency, and 
(cos a) the angular average of the scattering phase func- 
tion ( Bohren fc Huffmanl . 1983h . Q pr is dependent on the 
grain size. In Herbig AeBe stars with a radiation temper- 
ature of 10000 K, Q pl will be of order 1 for all grain sizes 
considered here, whereas in the cooler radiation field of T 
Tauri stars, it will be of order 1 for a > lpm, and some- 
what smaller for 0.1 pm grains. The grains will not move 
toward the star if the radiation force can at least balance 
the drag force of the gas streaming toward the star. The 
limiting case is reached if the grain is at rest with respect 
to the star. Then the gas streams by at a velocity v = Ug. 
As the free-fall flow is supersonic, we must apply the fric- 
tion law for the supersonic case in which we have 



To keep the grain fixed in space, the friction force must 
be balanced by the radiation force, 



F fric + F T 



ad 







(11) 



where we are neglecting the gravitational force on the dust 
grains, restricting the discussion to grains that are clearly 
inside the blow-out size regime. This equality is satisfied 
for a critical mass accretion rate 



M rr = 



4.8 x 10~ 8 



pr 



yr L Q 1 



Ms 



10000AU 



(12) 
(13) 



-Ffric = -Tta 2 pv 2 



(10) 



As we can see, the mass loss rate required to drag dust 
grains along is a function of distance from the star. Far 
away from the star, a relatively high mass infall rate will 
be required for the dust grains to be dragged along. As the 
infall moves closer to the star, it becomes easier to drag the 
dust grains along. Therefore, the true limiting infall rate 
is determined at large distances, where the free fall toward 
the star begins. If the infall starts at 10000 AU from the 
star, we find a minimum infall rate of 4.8 x 10~ 8 M Q /yr 
for solar parameters. 

It is therefore clear that gas drag may help to overcome 
the outward pushing radiation force only for the highest 
infall rates. For Herbig AeBe stars, the infall rates result- 
ing from our estimate contradict directly the observations 
of low Ay, and imply that low- level infall cannot play a 
role. For T Tauri stars, there is a window of infall rates be- 
tween approximately 4 x 10~ 7 and 5 x lO _8 M0/yr, which 
in principle would allow accretion of dust grains onto the 
disk. 

5.3. General remarks 

In this paper we have studied two scenarios to keep the 
disk optically thick in spite of dust coagulation. 

The first scenario is the low rate of matter infall onto 
the disk. This could help to retain a sufficiently high op- 
tical depth. However, it is unclear if these infall rates are 
in fact present. We can place severe limits on the allowed 
rates by considering the optical depth of the infalling enve- 
lope, and the effects of stellar winds and radiation pressure 
on dust (see Sect. I5.2[) . 

The second scenario involved beginning the coagula- 
tion simulation with a reduced abundance of fine dust. 
We have justified these special initial conditions with on- 
going coagulation during disk build-up. However, we have 
concluded that the special initial state needed to sustain 
T v ^ 1 without dust replenishment cannot be achieved in 
the time available. Even if this were achievable, the opti- 
cal depth that one can sustain for 1 Myr does not exceed 
Ty ~ 1, in contradiction with observational evidence. 

Therefore, fragmentation and erosion of aggregates ap- 
pear still to be the most promising explanation for the 
observed small dust grain populations in disks of 1 Myr 
age or older, as proposed by DD05. 
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It should be noted that the models we have presented 
in this paper are based on considerable simplifications. 
The most important of these is that we do not include the 
accretion of gas in the disk, toward the star. Any material 
that is located at 1 AU will be advected toward the star 
on a viscous timescale, which for a turbulence of a = 10 -2 
(assuming that the Schmidt number is 1) equals 10 4 years, 
and the structure of the disk may evolve significantly over 
a period of 1 Myr. As stated before, we have not included 
radial drift of grains. It is clear that the next step will 
have to be the full time-dependent modeling of the cou- 
pled equations of disk formation and evolution, and dust 
coagulation. With such models one could repeat the study 
completed here in a fully consistent way. 

6. Conclusion 

With this study, we arrive at the following conclusions: 

1. Low- level infall onto a disk can have a surprisingly 
strong influence on the steady-state abundance of 
small dust grains in a protoplanetary disk, in particu- 
lar if the infall reaches a disk that is initially cleared of 
small grains. This is due to a strong non-linear effect: 
the high initial dust content of the disk leads to very 
fast coagulation and a complete loss of small grains. 
The low dust densities produced by constant trickling 
are balanced by a far slower coagulation. 

2. To be significant, the trickling rate must exceed 
10~ n Mo/yr. Over a lifetime of lMyr, this corre- 
sponds to a total trickling mass in gas of 1O _5 M . 

3. In our calculations, the optical depth of the disk is 
produced by grains smaller than 3/im. In this case, we 
expect that the lO^im feature, an important indicator 
for small grains, will be present. Radiative transfer cal- 
culations will be required to determine whether this is 
the case. 

4. Almost perfectly laminar disks can remain marginally 
optically thick even without trickling, but for these 
disks the optical depth is produced by large grains (> 
3/im). 

5. Radiation pressure in Hcrbig stars excludes trickling 
as a viable mechanism for keeping the disk optically 
thick. Even for T Tauri stars, mass infall rates higher 
than 10~ 8 M Q /yr are needed to overcome the effects of 
radiation pressure on small grains. Large grains (above 
the blow-out limit) can be accreted more easily. For 
Herbig stars, the size limit of the order of a ^ 10/im, 
rendering it useless. In T Tauri stars, we find a ^ l^m, 
which is dependent on luminosity. 

6. Low initial dust-to-gas ratios could lead to disks with 
long-lasting optical depth of order unity. However, it 
appears unlikely that these states correspond to initial 
conditions because the time required to achieve these 
states is comparable to their duration. 

The above conclusions can be summarized in the fol- 
lowing way: While a low-level trickling infall of residual 



envelope matter onto the disk might conceivably sustain 
the disks 's optical depth above unity even without grain 
fragmentation, the conditions for this residual infall are 
very stringent. Ongoing and perpetual aggregate fragmen- 
tation and erosion in the disk as an explanation for the 
disk's optical depth remains a more effective explanation. 
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